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Limit to 2D mobility in modulation-doped GaAs quantum structures: How to achieve 

a mobility of 100 millions 

E. H. Hwang and S. Das Sarma 

Condensed Matter Theory Center, Department of Physics, 

University of Maryland, College Park, MD 20742-4111 

(Dated: April 29, 2008) 

Considering scattering by unintentional background charged impurities and by charged dopants 
in the modulation doping layer as well as by GaAs acoustic phonons, we theoretically consider the 
practical intrinsic (phonons) and extrinsic (background and dopants) limits to carrier mobility in 
modulation doped AlGaAs-GaAs 2D semiconductor structures. We find that reducing background 
impurity density to 10 cm~ along with a modulation doping separation of 1000 A or above will 
achieve a mobility of 100 x 10*^ cm^/Vs at a carrier density of 3 x lO" cm"^ for T = IK. At T = 4 
(10)K, however, the hard limit to the 2D mobility would be set by acoustic phonon scattering with 
the maximum intrinsic mobility being no higher than 22 (5) x 10^ cm^/Vs. Detailed numerical 
results are presented as a function of carrier density, modulation doping distance, and temperature 
to provide a quantitative guide to experimental efforts for achieving ultra-high 2D mobilities. 

PACS numbers: 73.40.-c, 73. 21. Ac, 71.30.+h 
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One of the most spectacular achievements of modern 
materials science is the continued enhancement of low 
temperature (T < IK) 2D electron mobilities in confined 
GaAs structures (e.g. quantum wells, heterostructures) 
over the last 30 yearsi^i^. In particular, the invention 
of modulation doping technique^ allowing spatial sepa- 
ration between dopants and carriers in 2D semiconduc- 
tor quantum structures and the continuous improvement 
in molecular beam epitaxy (MBE) eliminating uninten- 
tional charged impurities from the background have led 
to an exponential increase in the low-temperature 2D car- 
rier mobility in modulation doped GaAs-AlGaAs quan- 
tum structures from roughly 5 x 10'^ cm^/Vs in 1977^'^ to 
3.6 x 10^ cm^/Vs in 2007'', an amazing increase of mobil- 
ity by a factor of 7000 in 30 years. Although much less 
widely known than Moore's law, this is roughly equiva- 
lent to a factor of 3 increase in effective carrier velocity for 
every three years, more or less the same exponential rate 
of increase as in the celebrated Moore's law for Si micro- 
processors except that instead of miniaturization, which 
drives the Moore's law in Si electronics, modulation dop- 
ing directly enhances the low temperature carrier mo- 
bility. This astonishing 7000-fold low-temperature mo- 
bility enhancement in GaAs based 2D systems has been 
possible entirely through the systematic suppression of 
charged impurity scattering and interface roughness scat- 
tering in GaAs-AlGaAs heterostructures and AlGaAs- 
GaAs-AlGaAs quantum wells (which can be thought of 
as double heterostructures). This exponential increases 
in 2D mobility is therefore a triumph of materials science, 
not lithographic processing. The primary motivation for 
enhancing 2D mobility comes from the physics of frac- 
tional quantum Hall effect (FQHE) where higher mobil- 
ity leads invariably to the discovery of new phenomena*. 
Recent theoretical developments^ on topological quan- 
tum computation have provided added impetus^ to in- 
creasing 2D mobilities as high as possible. 

While improvement in materials science and device 



fabrication has continued to enhance the low tempera- 
ture (< IK) disorder-scattering limited 2D mobility, the 
question naturally arises on the possible intrinsic, rather 
than the extrinsic (i.e. limited by background or remote 
charged impurities, interface roughness, etc.) in quantum 
2D structures. This intrinsic limit, which must depend 
on the temperature, is obviously determined by phonon 
scattering since phonons, unlike impurities, can never be 
eliminated from a system except at T = OK. This in- 
trinsic temperature dependent mobility limit imposed by 
phonon scattering is sometimes referred to as the hard 
mobility limit in contrast to the soft limit imposed by im- 
purity scattering or interface roughness scattering or al- 
loy scattering which can, in principle, be reduced or even 
eliminated through materials science advances. The only 
way to reduce the phonon scattering is to reduce the tem- 
perature, and at any particular temperature, the phonon 
scattering indeed sets an impassable hard intrinsic limit 
to the mobility at that temperature. 

In addition to the obvious contributions by phonons 
to the intrinsic resistivity, there is actually an additional 
intrinsic contribution arising from the scattering by the 
charged dopants associated with the modulation doping 
producing the carriers themselves. This is a rather sub- 
tle effect that arises in semiconductors because an intrin- 
sic semiconductor does not have free charge carriers (as, 
for example, metals do with the metallic Fermi level ly- 
ing in the conduction band), and any mobile electron in 
a semiconductor must leave behind a positively charged 
dopant somewhere. The scattering between the dopant 
ions and the carriers must be considered an intrinsic resis- 
tive mechanism since a complete removal of the dopants 
would also eliminate the carriers themselves. Of course, 
the whole point or motivation for the modulation doping 
technique is to strongly suppress carrier scattering by 
the dopant ions through the large spatial separation be- 
tween the carriers in the 2DEG in GaAs and the dopants 
placed in a thin layer (so-called delta-layer) in AlGaAs 



at a distance d from the 2D carriers. Making the sim- 
ple assumption, which is accurate in this case, that the 
dopant ions and the 2D electron gas form together a par- 
allel plate capacitor, we get (e is the electron charge): 
n = K\AV\/{4:7T\e\d), where n, k, d, and \V\ are the 2D 
carrier density, the background lattice dielectric constant 
of the insulating regime (i.e. Gai-ajA^As), the modula- 
tion doping separation between the delta doping layer in 
AlGaAs and the 2DEG in GaAs, and the voltage or the 
energy barrier at the GaAs-Gai_2.Al2,As interface (which 
depends on the Al molar fraction a;) which creates the 
confining potential producing the 2DEG. We note that 
the maximum possible value of \AV\ « 0.8 eV for x = 0, 
and for a given x, the value of the modulation separation 
d determines the carrier density. In our discussion of 
intrinsic mobility we will ignore the precise relationship 
between d and n, and assume that d and n are indepen- 
dent variables as far as the 2D mobility calculation goes. 
This is mathematically justified because \AV\ could, in 
principle, be adjusted, making d and n independent vari- 
ables. 

The 2D conductivity a and mobility /x are given by 
a — nefi — ne^r /m where r is the scattering time 
or the transport relaxation time. There are many in- 
dependent contributions to the scattering rate t~^ in 
the modulation-doped GaAs- AlGaAs 2DEG system: (1) 
Acoustic phonon scattering via deformation coupling; 
(2) acoustic phonon scattering via piezoelectric cou- 
pling; (3) impurity scattering by unintentional back- 
ground charged impurities invariably present in the back- 
ground; (4) impurity scattering by intentional dopants in 
the modulation-doping delta layer; (5) interface rough- 
ness scattering at the GaAs-AlGaAs interfaces; (6) al- 
loy disorder scattering arising in Ga^Ali_a;As; (7) LO- 
phonon scattering via long-range polar Frohlich coupling; 
(8) short-range scattering by neutral defects and impuri- 
ties. Assuming Matthiessen's rule, which is strictly valid 
only at T = (but is reasonably valid for our purpose at 
higher temperatures also), we can write t~^ = Xiz ''"i" j 
where t^ indicates the transport relaxation time due to 
various individual scattering mechanisms. As discussed 
above, only the first two mechanism arising from the 
acoustic phonon scattering are true intrinsic resistive 
mechanism, but we will also include in our consideration 
the Coulomb scattering due to the intentional dopants 
in the modulation doping layer as well as the scattering 
by the unintentional random charged impurities invari- 
ably present in the background since in currently existing 
ultra-high mobility 2D samples, the background random 
charged impurity scattering is known to be the mobility- 
limiting mechanism. We note that the other scattering 
mechanisms (items 5 to 8 listed above, are known to be 
much less quantitatively important than the mechanisms 
(items 1 to 4) we are considering in this work. Polar 
carrier scattering by LO phonons is important in GaAs 
only at relatively high temperatures (> 50K), becoming 
dominant at room temperatures. Due to the rather high 
energy of the GaAs optical phonons (^ 36 meV), LO 



phonon scattering is completely suppressed in the tem- 
perature regime (< lOK) of interest to us in this work. 
Resistive scattering by optical phonons in 2D GaAs sys- 
tem has been considered in the literature^. 

We start by considering the Coulomb scattering by 
the 2D remote charged dopants in the modulation dop- 
ing layer which must invariably be present in modulation 
doped GaAs- AlGaAs structures and by the unintentional 
3D background charged impurities. Using the Boltzmann 
theory, the scattering rate, t~^, due to random charged 
Coulomb impurities, either in the background (i.e. 3D 
unintentional dopants) or in the modulation doping layer 
(i.e. in a 2D layer a distance d from the 2DEG), is given 

by8,9 
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where N^"'{z) is the 3D distribution of the a-th kind 

of charged impurities, u^ is the 2D Fourier transform 
of the 3D screened Coulomb impurity scattering matrix 
element for a 2D electron, scattering from the 2D wave 
vector k to k', 0kk' is the scattering angle between the 2D 
wave vectors k and k', and Sk = h?k^/2m is the 2D elec- 
tron energy. We use RPA to calculate screening in this 
work*'^. The two kinds of charged impurities important 
for our consideration are the 2D modulation dopants, 
Ni{z) = ndS{z—d), and the 3D unintentional background 
charge, Ni{z) = N^d, with rid and N^d denoting 2D and 
3D impurity densities respectively. To obtain the mo- 
bility, we obtain a thermal average of T~^(ek) over the 
finite temperature Fermi distribution function following 
the Boltzmann transport theory prescriptionSiS. 

For calculating the acoustic phonon scattering contri- 
bution to 2D carrier transport, we follow the Boltzmann 
theory developed by Kawamura and Das Sarmai^ii, con- 
sidering scattering through both the deformation poten- 
tial coupling and the piezoelectric coupling taking into 
account only the GaAs acoustic phonons. We do not 
show here the theoretical details, referring instead to the 
original workSiZ for the standard technical details. We 
include screening effects in our phonon scattering calcu- 
lations. 

Before presenting our calculated mobility results, we 
mention that our quantitative results depend on the pa- 
rameters d, Ud, N^D, and D, characterizing respectively 
the remote modulation dopant scattering {d, Ud), the 
background unintentional dopant scattering (N^o), and 
the GaAs deformation potential coupling constant (D) 
for the acoustic phonon scattering. One expects Ud = n, 
where n is the 2D carrier density, in the modulation 
doped structures, but for the sake of generality, we will 
assume n^ to be an independent variable, which is al- 
lowed in the presence of an external gate or compensating 
impurity charges. The precise value of the deformation 
potential coupling constant D to be used in 2D GaAs 
carrier transport consideration has been controversial as 
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FIG. 1: (Color online) (a) Calculated mobility times 2D impu- 
rity density (find) as a function of electron density for various 
values of d, and (b) exponents (a) of mobility in /i oc n" . 



discussed in refs. [y and|3, with values of D differing by a 
factor of two {D = 7 — 14eV) being quoted in the litera- 
ture. Since /i ^ |-D|^^, this is not a trivial matter as the 
calculated mobihty will be uncertain by a factor of four 
depending on the precise choice of D. We follow earlier 
transport theories, most notably^ and choose D — lOeV 
as the most suitable value for the GaAs acoustic phonon 
deformation coupling constant. Other values of D would 
reflect in an appropriate {D/10)^ rescaling of our calcu- 
lated phonon mobilities. 

All our results are for quantum wells of width 300 A 
which are symmetrically delta-doped on both sides with 
a modulation doping layer separation of d. We use 300 A 
wide symmetrically doped AlGaAs-GaAs-AlGaAs quan- 
tum wells for our calculations because recent ultra-high 
fabricated high-density 2D systems have mostly been 
fabricated in quantum well device (since higher carrier 
densities can be achieved in such symmetrically mod- 
ulation doped devices leading to higher carrier mobili- 
ties). All the qualitative trends we obtain remain valid 
for both modulation-doped AlGaAs-GaAs heterostruc- 
ture and quantum well systems with only minor quanti- 
tative variations arising from variations in system widths. 
Subband wavefunction effects^i^i are explicitly included 
in our theoretical calculations of various scattering ma- 
trix elements and the screening function, making our the- 
ory quantitatively reliable. We neglect all inter-subband 
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FIG. 2: (a) Calculated mobility times 3D impurity density 
(/i X Nzd) as a function of electron density. In inset the 
exponent a is shown. 



scattering, nonlinear screening, and localization effects. 

In Fig. 1(a) we show our T = mobility results for 
remote dopant scattering as a function of the electron 
density for various values of the modulation doping sep- 
aration distance d = — 1200 A (note that d = implies 
that the charged impurities are at the GaAsAs-GaAs in- 
terface). Since /i ex rf^ , we show in Fig. 1 the product 
[in^ as a function of n for various values of d. Our results 
at the lowest and the highest carrier densities are not re- 
liable due to our neglect of localization and intersubband 
scattering effects, respectively, but the results should be 
valid in the n = 10^° — 3 x 10^^ cm~^ range without any 
problem for high-quality samples. 

Assuming a high value of n^ = 10^^ cm~^, we see that 
the maximum intrinsic mobility at n = IQ^^ cm~^ in- 
creases from 3 X 10^ cm^/Vs for d = lOOA to the very 
large number of 30 x 10^ cm^/Vs for d = 1200A. Since 
the maximum achieved low-temperature 2D mobility so 
far has been around 36 x 10^ cm^/Vs for d > 1000 A, we 
conclude that the remote dopants are not the primary 
mobility limiting mechanism although they may be con- 
tributing about 30% of the resistive scattering. In Fig. 
1(b) we show that approximate density exponent (a) of 
mobility, which has often been discussed in the litera- 
ture, by writing fj, ^ n", i.e. a = dln/i/dlnn. We 
note that except for d = 0, a > 1 for remote impurity 
scattering, and a « 1.5 for larger values of d. This is 
simply a reflection of the fact that the remote scatter- 
ers are rather ineffectively screened by the 2D carriers, 
leading to a w 3/2 which corresponds to the unscreened 
limit. We conclude from Fig. 1 that, if all other scatter- 
ing mechanisms can be eliminated (most importantly, the 
unintentional background impurities), then at low tem- 
peratures, the 2D mobility could easily exceed 100 x 10® 
cmVVs at n = 3 X 10" cm"^ provided d > lOOOA. Note 
that for large d, where kpd ^ 1 with kp ~ -^n being the 
2D Fermi wave vector, /i ^ d"^, and therefore, as a matter 
of principle, remote dopant scattering can be reduced in- 
definitely simply by increasing d although in practice this 




FIG. 3: The mobility exponent a in /i cc n" for different 
Nch/Nba = 0, 1/5, 1/3, 1/2, 1, 2, 3, oo (from top to bottom). 
Here Nch (Nba) indicates the channel (barrier) 3D impurity 
density. The dashed line represents Nch = Ni,a- 



is difRcult since inducing carriers in the 2DEG becomes 
difficult for large d. 

In Fig. 2, we show our calculated mobility due to scat- 
tering by unintentional background charged impurities 
using the same quantum well parameters, but now with 
a three-dimensional uniform random background distri- 
bution of charged impurities of concentration N^o per 
unit volume which reside everywhere (both in GaAs and 
in AlGaAs). Again, fi ex iVg^, and therefore we show in 
Fig. 2 our calculated fiN^D as a function of the carrier 
density n at T = 0. For very clean and pure materials 
as typically used in the very best MBE growth labora- 
tories, N31J is extremely small, and cannot be directly 
measured by any spectroscopic tools (e.g. DLTS). As- 
suming N3D = 4 X 10^'^ cm^'^, which is a low estimate, 
we get ^ « 30 X 10*^ cm^/Vs at n = 3 x 10" cm'^. It 
therefore appears that the current limit to the achiev- 
able 2D mobilities in GaAs systems is set entirely by the 
background unintentional charged impurities, and reduc- 
ing their level to 10^^ cm~^ or below should raise the 2D 
mobility to 100 x 10^ cm^/Vs or above at low tempera- 
tures. To further reinforce the relevance of background 
scattering, we show as an inset to Fig. 2 the exponent 
a (i.e. fi ^ n") for background scattering by calculating 
a = d In /i/d In n. The value of a < 1 agrees well with the 
corresponding experimentally measured exponent. Note 
that the exponent a decreases as the width of quantum 
well decreases and the measured exponent in the ultra- 
high mobility system is considerably below a « 1.5 im- 
plied by remote impurity scattering [Fig. 1(a)]. 

Since the current low temperature {< IK) 2D mo- 
bility limit in high-quality GaAs structures is set by the 
Coulomb scattering by charged impurities, mostly by the 
unintentional background impurities and partially by the 
remote dopants in the modulation layer, we show in Figs. 
3 and 4 some more details on the theoretical mobility aris- 
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FIG. 4: (Color online) (a) Calculated mobility as a function 
of 5-layer distance d for fixed values of the densities n = 3 x 
lO^^cm"^, na = 10"cm"^, and N3D = lO^^cm"^ and (b) as 
a function of 3D impurity Nsd for n = 3 x 10 cm~ , Ud = 
10"cm"^, and d = 1200A. The dashed (dot-dashed) curves 
in (a) and (b) show remote (background) charged impurity 
contribution. 



ing from the screened Coulomb scattering. Noting that 
the quantitative details of the background unintentional 
doping are simply unknown (since their concentration is 
far too low to be studied spectroscopically), we show in 
Fig. 3, the calculated mobility exponent a (i.e. ^ '^ n") 
for different combinations of the barrier (i.e. AlGaAs) 
and the channel (i.e. GaAs) background impurity scat- 
tering leaving out the remote dopant scattering. The 
dashed curve with Nch/Nta = 1 in Fig. 3 indicates the 
results shown in Fig. 2 where the 3D random charged im- 
purities are distributed uniformly throughout the back- 
ground (i.e. both in the GaAs channel and in the Al- 
GaAs barrier). As expected, the exponent a decreases as 
the relative background scattering from the GaAs layer 
increases with respect to that from the AlGaAs barrier 
region. This is a direct effect of screening by the 2D car- 
riers themselves which is, of course, stronger (weaker) for 
the channel (barrier) impurities. For n « 10"'^"'^ cm~^, the 
exponent a ~ 0.7 is in excellent agreement with exper- 
imental results in the highest mobility variable density 
gated structures'^ where the 2D mobility has been mea- 
sured as a function of carrier density. It is noteworthy 
that in Fig. 3, a increases with increasing carrier density, 
reflecting the well-known^'^ peculiarity of 2D screening 
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that it becomes weaker (stronger) with increasing (de- 
creasing) carrier density since the relevant dimensionless 
screening parameter q^pjlkp^ where q^p (kp) are 2D 
Thomas-Fermi screening wave vector (Fermi wave vec- 
tor) goes as qTF/2kF ~ n~^/^ by virtue of qrp being 
density independent in 2D. Thus, for very large density, 
qTrf^kp <C 1, we recover the unscreened result a ^ 3/2. 
At high density {n :^ 10^^ cm~^), it is possible that in- 
terface roughness scattering (not included in our theory) 
could start contributing as the 2D electrons are pushed 
close to the GaAs-AlGaAs interface by the self-consistent 
confining potential, leading to a ^ since the interface 
roughness is a short-range scattering potential with very 
weak density dependence. 



In Fig. 4, we show the combined effects of remote and 
background Coulomb scattering (i.e. a combination of 
Figs. 1 and 2) on the 2D mobility at T = 0. In Fig. 
4(a), we show our calculated mobility as a function of 
the modulation separation d for fixed values of the den- 
sities n, Tid, and N^p). It is clear that for n^ = 10^^ cm^^ 
and N^p) — 10^^ cm^'^, /i reaches a maximum of around 
50 X 10^ cm^/Vs for d w 1200A since the background 
impurities become dominant for large d. Any further en- 
hancement of mobility would necessitate decreasing N^d , 
which is shown in Fig. 4(b), where we plot /i as a function 
of N^p) for fixed values of d, Ud, and n. A 2D mobility of 
100 X 10^ cm^/Vs (the '100 millions' in the title of our 
paper) can be achieved if N^d can be brought down to 
10"'^^ cm~^, an extremely small, but not an impossible, 
number. Results shown in Fig. 4 explicitly demonstrate 
how depending on the system parameters n, d, n^, and 
N^p, the actual mobility can be enhanced either by in- 
creasing d or by decreasing N^p , but extremely large 2D 
mobility would necessarily require large d (> lOOOA) and 
small N^p (< 10^'^ cm~'^). We show our intrinsic phonon 
scattering results next, emphasizing that for T < IK, the 
results of Figs. 1-4 shown above remain unaffected by 
phonon scattering. 

In Fig. 5 we consider the mobility due to acoustic 
phonon scattering, which is shown as a function of car- 
rier density for T — 4K and lOK. For lower temperatures, 
fiph increases by a large factor (/^ ^ T^^ for deformation 
potential scattering and fi ^ T^^ for piezoelectric scat- 
tering) since one is in the Bloch-Griiniesen range where 
phonon accupancy is suppressed exponentially^^. The 
important point of Fig. 5 is that at T = lOK (or even at 
4K), the mobility at n = 3 x 10^^ cm~^ is totally domi- 
nated by phonons and becomes 5 x 10^ cm^/Vs (22 x 10^ 
cm^/Vs), which is much lower than that for charged im- 
purity scattering. We therefore conclude that it will be 
impossible to raise 2D mobility above 22 x 10^ (5 x 10^) 
cm^/Vs at T = 10 (4)K since acoustic phonon scattering 
sets the intrinsic limit at these higher temperatures. At 
IK or below, however, acoustic phonon scattering is com- 
pletely suppressed, and the mobility is limited by back- 
ground charged impurity scattering and, to a lesser de- 
gree, by the remote dopant scattering. For example, at 
T — 2K the acoustic phonon limited mobility exceeds 
100 millions. 

Finally, as a conclusion, we combine the acoustic 
phonon and background impurity scattering in Fig. 6, by 
showing how progressive reduction of the unintentional 
background charge would raise the mobility toward 100 
million cm^/Vs provided d > lOOOA. We estimate that a 
background impurity density of 10^^ cm^'^ would produce 
a 2D mobility of 10^ cmVVs at T = li^ and n = 3 x 10" 
cm"2 foj. ^ ^ 1200A and n^ = 10" cm^^ 
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